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Moments  are  predicted  due  to  an  inviscid  liquid  payload  flowing  in  a  sequence  of  end  to 
uniform  cylinders  stacked  in  columns  displaced  off  the  spin  axis  of  a  coning  projectile  during 
free  flight.  These  moments  are  then  compared  to  moments  generated  by  flow  of  the  same 
liquid  that  saturates  porous  media  contained  in  a  sequence  of  uniform  cylinders  along  the 
projectile  symmetry  axis.  A  modification  to  the  classical  Stewartson  theory  and  is  added  to 
describe  inertial  waves  in  porous  media.  This  theory  is  used  to  analysis  the  inertial  waves 
produced  in  the  liquid  by  the  projectile  coning  motion.  The  resulting  side  moments  are 
examined  in  reference  to  the  type  of  porous  media  used  and  the  assumed  cylinder  geometry 
subjected  to  applied  coning  frequencies.  Eigen-frequencies  and  their  impact  on  liquid 
moments  are  discussed  concerning  the  stability  of  the  free  flight  projectile. 
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Nomenclature 


CT  Radial,  azimuthalporoucitycoefficent 

Cx  Axilporouscitycoefficient 

D  Cylinder  Diameter 

H  Cylinder  Length 

N  Number  of  sub-cylinders  per  candle-stick 

P  Projectile  angular  velocity  along  X'  axis 
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T  Non-dimensional  coning  frequency  TP  =  Frequency 
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Body  fixed  coordinate  system 

Sub-cylinder  aspect  ratio  A  =  2C/N 
Damping  rate  of  spiral  motion 
Fluid  mass  density 


Euler  roll,  pitch  and  yaw  angles  of  the  projectile 


I.  Introduction 

Predicting  the  moment  due  to  a  liquid  payload  in  a  spinning  and  coning  projectile  is  a  problem  of  considerable 
interest  to  Army.  Stewartson1  considered  the  linear  problem  of  a  liquid  payload  in  a  spinning  right  circular  cylinder 
using  separation  of  variables  and  eighenvalue  expansions  for  an  inviscid  liquid,  i  First  order  viscous  boundary  layer 
corrections  to  the  Stewartson1  theory  were  carried  out  by  Wedemeyer2  and  Murphy.3  A  method  for  calculating  the 
linear  liquid  moment  using  the  full  linear  viscous  equations  with  boundary  layer  corrections  confined  only  to  the 
endcanps  was  also  presented  by  Flail  and  Sedney  and  Gerber4'5. 

A  further  interest  to  the  army  is  to  consider  a  series  of  uniform  circular  cylinders  stacked  end  to  end  separated  by 
impenetrable  end-caps  (candle-sticks).  These  candle-stick(s)  may  be  situated  along  the  symmetry  axis  or  off  set  but 
parallel  to  the  symmetry  axis  of  the  projectile.  Coning  motion  induced  liquid  moments  are  considered  here  for  a 
number  of  candle-stick(s)  configurations.  The  eigen-frequencies  for  such  configurations  are  shown  to  be  identical  to 
those  found  by  Stewartson1. 

Liquid  payloads  contained  in  a  highly  permeable  material  have  also  been  of  interest  to  the  Army  for  some  time. 
Laboratory  tests  and  flight  tests  have  shown  that  a  highly  preamble  medium  can  significantly  reduce  the  spin-up 
time  of  a  liquid  payload.6  Flight  stability  for  liquid  saturated  permeable  payloads  has  also  been  examined  by 
D’Amico.7  The  work  here  extends  the  Stewartson1  problem  by  considering  a  cylindrical  cavity  filled  with  a 
permeable  medium  that  is  impregnated  with  an  invisicd  liquid.  A  further  modification  is  introduced  by  segmenting 
the  cavity,  along  the  symmetry  axes,  into  a  sequence  of  equal  length  cylinders.  Each  of  these  cylinders  is  separated 
by  impermeable  end-caps.  The  porous  media  is  modeled  by  a  drag  term,  which  is  proportional  to  the  liquid  velocity 
relative  to  the  assumed  ridge  porous  media  that  is  added  to  the  linearized  Euler  equations.  This  analysis  examines 
the  induced  liquid  moment  as  a  function  of  parameters  found  by  Stewartson  1  plus  parameters  describing  the  porous 
media  and  the  number  of  segments  in  the  cylindrical  cavity. 


II.  Equations  of  Motion  for  the  Candle-Stick  Configurations 


Figure  1.  shows  the  X',Y',Z'  axes 


Figure  1  Coordinate  Systems  of 
configuration 


rotating  uniformly  about  Z'  with  angular  speed  P  =  (P,  0,  0)  .  The  liquid 
is  assumed  initially  to  be  rotating  as  a  rigid  body  with  the  same  angular  speed 
P  so  the  velocity  Y'  of  the  liquid  inside  one  the  cylinders  is 

V'  =  Px(x,  R0cosB  +  rcos9,  R0sinB  +  rsin9)  (1) 

The  unperturbed  state  for  Eq.  (1)  satisfies  the  Euler  equation: 

p 

P2exx(exx(x,  R0  cosB  +  rcos0,  R0 sinB  +  rsin9))  = -V  — 

P 

,  (2) 

for  which  Ps  is  the  unperturbed  liquid  pressure.  Following  Stewartson  and 
letting  R  =  (x,  R0  COS  B  +  r  COS  9,  R0  sin  B  +  r  sin  0)  gives  the  scalar 
potential 
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(3) 


(4) 

(5) 


(6) 


Ps/p  =P2  R0  cos  (B  -  0)-  P2  r2/2 

Now  consider  a  perturbed  angular  velocity  of  the  projectile  for  small  coy,coz  given  by 

Q=P  +  (0,  C0y,  C0Z) 

and  adopting  projectile-fixed  X,  Y,  Z  axes  causes  the  liquid  velocity  to  take  the  form 

V  =  Px(x,  R0cosB  +  rcos9,  RosinB  +  rsin0)  + v 

where  the  components  of  V  have  the  same  order  of  magnitude  as  C0y  ,  C0Z  . 

The  Euler  equations  for  small  perturbations  written  in  the  body  fixed  frame  now  read  as 

dt  p  dr 

— +  2PV  +  — M  =  o 
dt  pr  d0 

dU+idp=0 

dt  p  dx 

where  (V,  W,  U)  are  the  cylindrical  components  of  the  perturbed  velocity,  p  is  the  perturbed  pressure. 

The  boundary  conditions  at  the  solid  wall  satisfy 

(v,  W,  U)*n  =  fixRs»n  (7) 

where  n  is  a  outward  unit  vector  on  the  wall  and  Rs  a  point  on  a  cylinders  wall.  Following  Murphy/Cooper  the 
perturbation  quantities  coy ,  wz  are  assumed  to  represent  coning  motion  given  by 

CO  =  K0  P  (s  T  sin(P  T  t)  +  (T  - 1)  cos(P  T  t))  eE  P  T  ‘ 

(8) 

coz  =  -K0  P  ((T  -  l)sin(P  T  t)  -  s  T  cos(P  T  t))esPTt 

where  the  coning  damping  rate  is  8  the  coning  frequency  is  T  and  K0  is  the  magnitude  of  the  coning  angle. 
Equation  (8)  is  now  used  in  Eq.  (7)  so  the  boundary  conditions  now  become: 

V  =  jRe(x  K0  P  (S  -  i)  ePSt+i9) 

U  =  jHe(  K0  P  (eiB  R0  +  r  ei6)(S  -  i)  ePSt+ie) .  (9) 


i  =  -y- 1’  S  =  (£  +  i)T 

These  boundary  conditions  suggest  Eq.  (6)  is  separable  giving 


"vl 

V 

w 

w 

u 

u 

_p_ 

_p_ 

PSt 


(10) 


and  solving  for  the  velocity  components  yields: 
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—  rS  +  2  — 


V  = 


dr 


d9 


w  = 


rpp(s2  +4) 

—  S-2r— 
d0  dr 

r  p  P  (s2  +  4) 
dp 


(11) 


u  = 


dx 

pPS 


Using  the  continuity  equation,  V  (v,  W,  u)  =  0  ,  produces  the  following  equation  for  the  pressure  p 


.^-P.UPU-P  -2_2dP=0 


-  +  r- 


-r  cr 


dr2 


dr  d0 


dx 


2  S  +  4 
a  = - 


(12) 


At  this  point  in  the  analysis  it  is  useful  to  consider  each  cylindrical  candle-stick  ranging  from,  —  C  <  X  <  C ,  to 
consist  of  a  end  to  end  sequence  of  N  equal  length.  A,  cylinders  with  impentetrible  end  caps  such  that 
A  =  2  C/N .  Applying  Eq.  (9)  to  each  of  the  sub-cylinders  shows  that  separation  of  variables  gives 

p  =  Aj  K0  P2  a2  p  cos  (n  (2j  +  1))(C  +  x)/A)  Jj  (n  (2j  +  l))a  r/A)e10 
+  K0  P2a2p(rxE/a2  +  rD/a)eI0  +  KoP2apFx 


(13) 


where 


8(-l)nA2(S-i)Z(S  +  2i) 


o2, 


A.  = 


j  1  (rr  (2  j  +  1 )  aa  /A)a(S  -  2i) 

n  a  (2j  +  l)S  a  jQ(7i  (2j  +  l)a  a/A) 
E  =  (S-i)S 

D  =  ((2n  +  l)A-2C)(S-i)2 
F  =  elB  R0(S-i)s/a 


2„ 2^,  ,  a2 


n  a 


(2j  +  lf 


(14) 


0  <  n  <  N,  j  =  0,1,2  • 


III.  Candle-Stick(s)  Liquid  Moments 

The  moment  induced  by  the  liquid  contained  in  the  segmented  cavity  is  calculated  from  the  time  derivative  of 
the  angular  momentum  field.  Non-dimensionalizing  the  moment  with  2 71  pa4  CP"  it  is  convenient  to  write  the 
side  moment  components  as3 

Mv  +iMz  =  xCLM  (2jipail  CP2)K,ePTI 
CL„=CLSM(T,N,C)+iCL1M(T.N,C) 

Unit  vectors  of  the  body-fixed  cylindrical  coordinates,  (ex,  er,  ee)  are  written  as  the  complex  quantities  in 
terms  of  (ex,  eY,  ez) 
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er  =  eYcos0  +  ezsin9  <=>  exp(i0) 

/  \  (16) 
ee  =  _eYsin0  +  ezcos0  <=>  iexp(i0) 

Placing  these  in  the  moment  integral  and  using  the  Reynolds  Transport  Theorm10  yields  the  following  expression  for 
the  liquid  moment  coefficient: 

T Clm  =  ~t p  X  f  f  fxex  +  r  er)  x  [(S  -  i)q  +  2  ex  x  q]-  i[r2  -  2  x]  ex jrdrdx 

2ji-n=oJ/  (17) 

a 

q  =  (v,  w,  u) 

The  last  result  written  terms  of  Eq.(l  1)  is  : 


TC,„=P^(S-i)J(S  +  2i) 


jt'N2 


^S-2i  2 

- -  —  +  1 

v  ct2S  ct2  j 


jt(zz)/jo(z) 


(ji(2j  +  l)aNS)- 


2  C(S  -  2i)j,(zz) 

Jo(zz) 


1  (S  -  i)((4  C2  +  3)N2-16C2)s2+16C2(l-N2)- 

2  i  (l  6  C2  -  N2  (3  +  8  C2)) 

12  N2  (S  —  2  i) 

+  iR2  S(S  -  i)sin B elB  ; zz  s  7t(2j  +  1)ctN/2C 


(14) 


(15) 


+  ie  lB (2ssinB  +  cosB)Rq  +- 


For  a  not  a  resonance  value  and  as  the  number  of  cylinders  becomes  large  (N— >co)  the  value  of  CLM 
approaches  the  frozen  liquid  limit  given  by: 

CLM  -»  -i e~iB  (2 s sin B  +  cos B)R2  T  -  (4 C  ~ 8  +  ^ T 

(4  C2  +  l)(s  +  i) 

4 

The  eigen-values  for  this  problem  can  be  found  by  inspecting  Eq.  (14)  and  after  some  algebra  these  are  found  to  be 
zeros  of 

Jj(zz)  ±  V  B  +  1-1 

zzJ0(zz)  B 

B  =  4zz2C2 
7t2(2j  +  l)2N2 

which  gives  the  Stewartson1  values  when  zz  is  solved  for  frequency  T  . 


(16) 


IV.  Equations  of  Motion  for  the  Axial  Porous  Media  Configuration 


For  this  problem  the  moment  arm  R0  =0  in  Eqs.  (3  &  5)  and  the  position  vector  becomes 
R  =  (x,  rcosO,  r  sin  0) since  the  candle-stick  axis  is  the  symmetry  axis  of  the  projectile.  However,  the  Euler 
equations  are  now  modified  to  account  for  an  inviscid  liquid  flowing  through  porous  media.  The  modification  is 
assumed  to  have  terms  proportional  to  the  liquids  velocity  relative  to  that  of  the  porous  media  which  is  taken  to  have 
the  same  velocity  as  the  coning  projectile.  The  Euler  equations  are  now  written  as: 
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dV 


2PW  +  C  p(v-v)+-—  =  0 

dt  p  dr 


dW 

dt 


+2PV+C 


dU 


t  p(w  -  w)+ 


1  dp 
pr  d0 


—  ^  =  0 


+  cxp(u-u)+-^-  =  o 

dt  p  dx 


where  the  velocity  of  the  media  obtained  from  the  rotation  kinematics  is  given  by 

V  =  xK0  P  (S  -  i)e  PSt+ie 
W  =  ixKOP(S-i)ePSt+i0 
U  =  -rK0P(S-i)ePSt+ie 

Separating  variables  according  to 


“v“ 

V 

w 

w 

u 

u 

_p_ 

_p_ 

pst+ie 


leads  to  the  solution  of  Eq.  (16) 


2ip  +  — (S  +  Ct)  , 

dr  xCtP(S-i) 

V  = - 1 - vH - - — ^ 

rpP((S  +  Ct)2  +  4j  S  +  Ct-2i 

i(S  +  Ct)p-2  — 

dr  ixCtP(S-i) 

rpP((S  +  Ct)2+4)  S  +  Ct-2i 

—  d^  +  rCxP(S-i) 
pP  dx 

u  =  - 


(17) 


(18) 


(19) 


(20) 


s  +  cx 

Continuing  with  the  continuity  equation  gives  the  following  equation  for  the  perturbation  pressure  p 

2d2p  dp  2_2d2p  A 

r  — -  +  r  —  -p-r  a  — -  =  0 

dr2  dr  dx2 


—2 

a  =  — 


(S  +  Ct)2  +  4 


(21) 


(s+c.Xs+cJ 

Solving  and  using  the  demands  of  Eq.  (17)  produces  the  following  coefficients  for  expression  Eq.  (13) 
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A;  = 


8  (-  l)n  A2  (S  -  i)2  (S  +  Ct  +  2i) 


"jI(7i(2j  +  l)aa/A)A(S  +  Ct-2i)- 

n  (2j  +  1)(S  +  Ct)o  j0(7t(2j  +  l)aa/A)y 
E  =  (S-i)S 

D  =  -((2n  + 1) A  -  2  C)(S  -  i)2 
F  =  0 


n 


(2j  +  l  f 


(22) 


0  <  n  <  N,  j  =  0,1,2  •  •  • 

V.  Porous  Media  Liquid  Moments 

Using  the  above  procedure  for  calculating  liquid  side  moments  results  in  the  following  expression  for  the  porous 
media  moments 

TCl„  =  -  !28l(:’  -  — (s - i)!(S  +  Ct  +  2 i) 


7t4a2N2(S  +  CtXs  +  Cx) 

2S3  +  2(CX  +  3Ct  -  i)S2  +  2(2CtCx  +  3C2  -  2iCt  +  4)s 
2(C,+C,-i)(c=+4) 

^2Jj(zz)(S  +  Ct)(S  +  Ct  -  4i)C  - 


XJiM/(2J  +  I) 


7x(2 j  +  l)aN  J0(zz)(S  +  CjS  +  Ct  -  2i) 


+ 


iC2s((N2 -2^-i(3N2  -4))/3N2- 


3N2(S  +  Ct  -  2i) 


+ 


iS(S  -  i)/4  -  2iC2(2Cx  +  Ct)/3N2(S  -  2i); 
zz  s  tt:(2 j  +  l)aN/2C 

Here  again  the  limiting  case  with  N  — >  oo  gives  the  value  of  CLM  for  the  frozen  liquid  side  moment 

TClm  ^  is((4C2  +  3)s  -  3i(4C2  +  l))/l2 

The  eigen-values  for  this  problem  are  found  from  Eq.  (22)  and  after  some  algebra  they  are  zeros  of 


J,(zz)  _±(Cx-Ct  +  2i)VB((Cx-Ct)2B-16)-16  |  (Cx-Ct)B-4i 
zzJ0(zz)  2((CX  -  Ct)(Cx  -  Ct  +  4i)-  16)B  2(Cx  -  Ct  +  2i)B 


(23) 


(24) 


(25) 


This  equation  shows  when  Cx  =  Ct  yields  eigen-values  ZZ  are  the  same  as  those  found  by  Stewartson',Eq.  (16), 
provided  zz  =  zz  . 
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VI.  Calculation  Method 


The  equations  of  the  last  sections  need  to  be  calculated  for  a  wide  range  of  flight  and  porous  media  parameters 
all  of  which  need  the  values  of  Bessel  functions.  For  small  values  of  |zz|  simply  using  power  series  expansions  of 

each  Bessel  function  works  very  well.  Bessel  functions  at  large  values  of  |zz|  were  obtained  by  using  asymptotic 
expansions  for  each  Bessel  funcion9.  Generally  calculating  Bessel  functions  for  complex  arguments  for  intermediate 
values  of  |zz|  is  a  non-trivial  problem  and  the  methods  used  here  employs  Gauss  continued  fractions.  This  author 

has  judged  that  a  further  discussion  of  these  methods  is  not  appropriate  for  this  article  but  the  reader  should  be  aware 
of  the  numerical  difficulties  associated  with  calculating  Bessel  functions. 


VII.  Results 

Figures  2  and  3  show  the  side  moment  CLSM  and  in  plane  moment  CLIM  as  functions  of  the  non- 
dimensional  frequency  T .  These  moments  show  the  eigen-frequencies  for  a  candle  stick  that  is  displaced 
off  the  symmetry  axis  by  R0  =1.5. 


Figure  2  Liquid  Side  Moment 


In  Plane  Moment  N=1 


Figure  3  Liquid  In-Plane  Moment 
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The  next  two  plots  given  in  Fig.  4  &  5  present  the  liquid  moments  for  three  candle  sticks  uniformly  distributed 
around  the  projectile  symmetry  axis  forR0  =  1.5 .  The  value  of  8  and  N  are  chosen  so  that  Stewartson1  eigen¬ 
values  are  not  present  for  the  given  range  of  T  .  Figure  5  also  displays  the  frozen  liquid  moment  which  indicates 
that  the  liquid  behaves  like  frozen  liquid  for  increasing  values  on  N  provided  eigen-frequencies  are  avoided. 


Figure  4  Sum  Of  Three  Side  Plane  Moments  Figure  5  Sum  of  Three  In  Plane  Moments 

Next  consider  a  centrally  located  candle  stick  containing  a  liquid  in  porous  media.  Figures  6  and  7  presents 
examples  of  liquid  moments  for  typical  values  of  CT  =  0.3,,  Cx  =  0.3, 0.5  and  8  =  —0.612,—  1.807  .  Once 

again  the  moments  display  the  eigen-frequency  behavior  of  saturated  porous  media  with  increasing  N  =  1 .  These 
indicate  that  porous  media  forces  the  eigen-frequencies  to  be  complex  when  ever  8  is  chosen  to  cause  <7  to  be  an 
eigen-value 


Side  Moment  With  Resonance 
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Figure  6  Side  Moment  For  Central  Porous 
Cavities 


In  Plane  Moment  With  Resonance 
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Figure  7  In-Plane  Moments  For  Central  Porous 
Cavities 


VIII.  Conclusions 

The  off  axis  candle  stick  problem  has  been  shown  to  be  equivalent  to  the  inviscid  Stewartson1  problem  whenever 
porous  media  is  not  present  or  can  be  ignored.  Resonances  are  independent  of  the  candle  stick  off  axis  position. 

In  the  case  of  a  candle  stick,  containing  saturated  porous  media,  located  along  the  symmetry  axis  of  the  projectile 
generally  forces  resonances  only  for  complex  values  of  S  .  In  all  cases  the  liquid  moments  approach  the  values  for  a 
frozen  liquid  when  eigen-frequencies  are  not  present.  In  the  particular  case  where  CT  =  Cx  it  is  possible  find  the 
resonances  from  the  Stewartson1  tables  but  if  CT  ^  Cx  requires  a  numerical  search  for  resonances  in  the  complex 
plane  for  ZZ  from  which  S  can  then  be  calculated. 
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